Astronomy & Astrophysics manuscript no. Marti- Vidal-optimize-LOFAR 
May 5, 2010 



©ESO2010 



Optimum estimate of delays and dispersive effects in 
low-frequency interferometric observations 



I. Marti- Vidal 



o 



J2 
Or 

6 

S-H 
•«— > 

CO 



> 

in 

m 

q 

in 

o 
o 



Max-Planck-Institut fiir Radiastronomie, Auf demHiigel 69, D-53121 Bonn (Germany) e-mail: imartiv@mpifr.de 

Submitted on 2009 Dec 23. Accepted on 2010 May 3. 

ABSTRACT 

Context. Modern radio interferometers sensitive to low frequencies will make use of wide-band detectors (with bandwidths of the 
order of the observing frequency) and correlators with high data processing rates. It will be possible to simultaneously correlate data 
from many sub-bands spread through the whole bandwidth of the detectors. For such wide bandwidths, dispersive effects from the 
atmosphere introduce variations in the fringe delay which change through the whole band of the receivers. These undesired dispersive 
effects must be estimated and calibrated with the highest precision. 

Aims. We studied the achievable precision in the estimate of the ionospheric dispersion and the dynamic range of the correlated 
fringes for different distributions of sub-bands in low-frequency and wide-band interferometric observations. Our study is focused on 
the case of sub-bands with a bandwidth much narrower than that of the total covered spectrum (case of LOFAR). 
Methods. We computed the formal statistical uncertainty of the ionospheric delay, the delay ambiguity and the dynamic range of the 
correlated fringes using four different kinds of distributions of the sub-bands: constant spacing between sub-bands, random spacings, 
spacings based on a power-law distribution, and spacings based on Golomb rulers (sets of integers, n { , whose sets of differences, 
rij - lit, have non-repeated elements). 

Results. We compare the formal uncertainties in the estimate of ionospheric effects in the data, the ambiguity of the delays, and the 
dynamic range of the correlated fringes for the four different kinds of sub-band distributions. 

Conclusions. For a large number of sub-bands (> 20, depending on the delay window) spacings based on Golomb rulers give the most 
precise estimates of dispersive effects and the highest dynamic ranges of the fringes. Spacings based on the power-law distribution 
give similar (but slightly worse) results, although the results are better than those from the Golomb rulers for a smaller numbers of 
sub-bands. Random distributions of sub-bands result in relatively large dynamic ranges of the fringes, but the estimate of dispersive 
effects through the band is worse. A constant spacing of sub-bands results in very bad dynamic ranges of the fringes, but the estimates 
of dispersive effects have a precision similar to that obtained with the power-law distribution. Combining all the results, the power- 
law distribution gives the best compromise between homogeneity in the sampling of the bandwidth, precision in the estimate of the 
ionospheric dispersive effects, dynamic range of the correlated fringes, and ambiguity of the group delay. 

Key words, atmospheric effects - techniques: interferometric - instrumentation: interferometers - telescopes 






1. Introduction 

Modern radio interferometers sensitive to low frequencies, like 
the LOw Frequency ARray (LOFAR, see, e.g., de Vos, Gunst, & 
Nijboer 2009 and references therein), will make use of wide- 
band detectors (with a bandwidth of the order of the observ- 
ing frequency) and correlators with a high data processing rate. 
For the case of LOFAR, it will be possible to correlate a band- 
width of 42 MHz, which can be divided into many sub-bands 
to be spread through a wide band (either 10-90 MHz, 1 10— 
190 MH z, 170 -230 MHz, or 2 10-250 MHz; de Vos, Gunst, & 
Nijboer 2009). Therefore, the correlator will be able to si- 
multaneously cover a large portion of the low-frequency spec- 
trum. However, atmospheric dispersion may introduce strong 
frequency-dependent effects, which will be so different for each 
sub-band and must be estimated and calibrated with care. At 
low frequencies, the ionosphere is the main limiting factor in 
the quality of the interferometric observations. For the calibra- 
tion of the ionospheric dispersion, the Total Electron Content 
(TEC) must be estimated through the field of view (FoV) over 
each antenna and for each time. The TEC can be estimated, for 
instance, from Global-Positioning-System (GPS) data (e.g. Ros 
et al. 120001 Todorova, Hobiger, & Schuh l200"6l With these tech- 
niques, however, it is not possible to obtain high-resolution esti- 



mates of the ionospheric turbulent screen over the FoV. For that 
purpose, it is necessary to apply self-calibration strategies, using 
therefore the interferometric data to derive the structure of the 
TEC screen over each station and its evolution (see Interna et 
al. [2009 and Cohen & Rottgering 2009 for a discussion on the 
ionospheric screens). 

The distribution of sub-bands through the whole bandwidth 
of the detectors affects the scientific information that can be ex- 
tracted from the observations, but is also important for a precise 
and accurate estimate of the atmospheric dispersion from the in- 
terferometric observables. The distribution of sub-bands has also 
an effect in delay space, since the interferometric fringes have a 
shape related to the Fourier transform of the bandpass (i.e., in 
our case, the distribution of sub-bands). Therefore, we should 
distribute the sub-bands in such a way that the spectral cover- 
age and sampling are maximized, but the dynamic range of the 
fringes (i.e., the height of the fringe peak relative to that of the 
highest sidelobe) is also maximized to improve the sensitivity 
of the interferometer. Finding out the right distribution of sub- 
bands to achieve an optimum spectral sampling and fringe dy- 
namic range is not that simple, and the answer may depend on 
each case, namely, the total bandwidth to be covered, the number 
of sub-bands, and the bandwidth of each sub-band. For instance, 
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Petrachenko (2008) studied the performance of "broadband de- 
lays", which are computed from several bands (up to 5) spread 
from 2-3 GHz to ll-14GHz, as a function of the way these 
bands are distributed through the spectrum. Petrachenko (2008) 
concluded that the use of more than 2 bands, covering a total 
bandwidth as wide as possible, improves the performance of the 
interferometer. 

In this paper, we report on a study of the spectral coverage, 
the precision in the estimate of atmospheric dispersive effects, 
the delay ambiguity, the dynamic range of the fringes, and the 
precision in the estimate of the source spectral index for differ- 
ent kinds of sub-band distributions and spectral configurations 
of an interferometer at low frequencies. The remainder of this 
paper is structured as follows: in Sect.|2]we describe the process 
of analysis. In Sect. [3] we describe the different sub-band distri- 
butions studied. In Sect. |4]we report on the results obtained and 
in Sect. |5] we summarize our conclusions. 



2. Analysis 

2.1. Spectral configuration of the interferometer 

The spectral configuration of an interferometer can be charac- 
terized using the following parameters: i) minimum observing 
frequency, v m , ii) total bandwidth in units of the minimum ob- 
serving frequency, i.e., 



P - Om - v m )/v m , 



(1) 



where vm is the maximum observing frequency, iii) number of 
sub-bands, N (and their distribution), and iv) the bandwidth of 
each sub-band, Av. The /th sub-band is centered at frequency v, 
(with v m < Vj < vm), which can be written as 



V,- = (l+j8*i)Vm, 



(2) 



being /?, a real number between Q.5Av/(J3v m ) and 1 - 
0.5Av/(j6v m ). The spectral configuration of an interferometric 
observation is then characterized by v m , [5, Av,, and /?,-. For sim- 
plicity, we assume the same Av, = Av for all sub-bands, and 
we also assume that Av << j3v m (i.e., the bandwidth of the sub- 
bands is much narrower than the total covered bandwidth, so /?, 
is defined between and 1). This latter assumption corresponds 
to the case of the LOFAR interferometer. 



2.2. Dynamic range of the fringes 

The response of a baseline of an interferometer to a set of ob- 
served sources is equal to the addition of several fringes, one 
fringe for each source. The amplitude of these fringes is equal to 
a shrinked version of the amplitude of the Fourier transform of 
the bandpass (e.g., Thomson, Moran & Swenson 1991). Spectral 
effects in the sources are not considered here. In the case of sev- 
eral sub-bands spread through a wider band, different shapes can 
be obtained for the fringes, which may have relatively large side- 
lobes. These large sidelobes may lead to confusion in the esti- 
mate of the fringe peaks. In frequency space, this effect can be 
understood as the possibility of fitting different slopes of phase 
vs. frequency to the same dataset. Minimizing the height of the 
sidelobes and maximizing their distance to the fringe peak de- 
creases the probability of confusion and, therefore, enhances the 
sensitivity of the interferometer. 



For each of the studied band distributions (described in Sect. 
|3}, the shape of the fringes, F, was computed as a vector with 
elements given by 



F k = 



N , 

> cos \n — — (k - 1) 



(3) 



i.e., the module of the Direct Fourier Transform (DFT) of the 
sub-band distribution (the sine term of the DFT is zero). The 
vector F is thus computed assuming that the bandwidth of the 
sub-bands, Av is much narrower than the total bandwidth, /3v m . 
The dynamic range of the fringe is estimated as the ratio between 
the fringe peak (which corresponds to the element Fj) and the 
peak of the highest sidelobe. We call this ratio D. The ambigu- 
ity of the delay is computed as the distance between these two 
peaks, in units of the Nyquist time resolution (i.e., 1 /(2ySv m )). We 
call this quantity r am b. For the results reported in the following 
sections, we used a vector F of 512 elements. 

If the width of the sub-bands is not much narrower than the 
total covered bandwidth, the dynamic range, D, is corrected by 



/ Av \ 

D -> D/Sinc 7T— — Tamb , 

\ 2j3v m ) 
where sinc(^) = sin (x)/x. 



(4) 



2.3. Estimate of the atmospheric dispersion 

The ionosphere introduces a change in the phase of the visibili- 
ties. For a given baseline and source, this phase depends on fre- 
quency as (e.g., Thomson, Moran & Swenson ll99H 



A0, = .Kv~ 



(5) 



The parameter K is related to the TEC of the ionosphere over 
the elements of the baseline in the line-of-sight direction to the 
source. For a good calibration of the ionospheric dispersion, K 
must be precisely estimated for each baseline, time, and pointing 
direction over the Fo V of the interferometer. 

The formal statistical uncertainty in the estimate of K is that 



of the slope of the linear fit of < 
show that this uncertainty is 



vs. v . It is straightforward to 



cr(K). 



(< y; 2 > 



< vj 1 > 2 r' 



crCv" 1 )- 1 , 



(6) 



i.e., the precision in K is proportional to the standard deviation 
of the distribution of the inverse of the central frequencies of 
the sub-bands, <r(v~ l ). This standard deviation maximizes when 
half of the sub-bands gather close to v m and the other half gather 
close to vm- In this case 



0"(V; ^max = 0.5 (vj - Vj^). 



(7) 



However, this distribution of sub-bands is not, by far, op- 
timum, since the sampling of the total band is very poor and, 
moreover, spectral effects in the sources, which would not be 
well sampled through the bandwidth, could introduce important 
systematics in the estimate of the ionospheric dispersion using 
Eq. Additionally, there could be several undetermined 27r- 
cycles of the phase drift caused by the ionosphere between v m 
and vm, so, for this distribution, it would not be possible to con- 
nect the phases between all sub-bands to obtain a correct esti- 
mate of K. 

The uncertainties in the estimate of the ionospheric disper- 
sion, which we analyze in Sect. [4] were estimated as cr(K) 
(computed from Eq. |6]l in units of its minimum possible value, 
o~(K) m in (i.e., that one corresponding to cr(yj l ) max , which is 
given in Eq.[7J. 
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Fig. 1. Examples of the four kinds of distributions studied in this paper (using 16 sub-bands and setting f3 - 0.5, see Eq. Q}. (a) 
corresponds to the uniform (i.e. constant) distribution, (b) to the random distribution, (c) to the power-law distribution, and (d) to 
the distribution based on the Golomb ruler. 



3. Distributions of sub-bands 

3.1. Uniform (i.e., constant) distribution 

This is the most simple spectral configuration of the interferom- 
eter. The central frequencies of the sub-bands are distributed as 



1+J3 



i-\ 

N -1 



(8) 



i.e., the frequency spacing between sub-bands is constant. We 
show an example of this distribution in Fig.Qla). 



3.2. Random distribution 

In this case, the sub-bands are randomly distributed over the 
bandwidth, i.e, 



(l+0Ui)v n 



(9) 



where £/, is a random real number between and 1, following 
a uniform statistical distribution. We show an example of this 
distribution of sub-bands in Fig. Qlb). Different sets of £/,■ may 
translate into different fringe dynamic ranges and precisions in 
the estimate of ionospheric dispersion (usually, a higher preci- 
sion in the estimate of the atmospheric dispersion translates into 
a lower dynamic range of the fringes). We estimated the quanti- 
ties cr(K), D, and r am b as the averages of those computed from 
100 different sets of [/,-. Therefore, the results reported in Sect. [4] 
for the random distribution correspond to an averaged behavior 
of this distribution. 



3.3. Power-law distribution 

The ionosphere introduces larger phase drifts at lower frequen- 
cies. Therefore, it is plausible that a distribution that samples 
better the region of lower frequencies will give more precise es- 
timates of the ionospheric dispersion, since the phase drifts will 
be better sampled in the region of the spectrum where the iono- 
spheric effects are larger. A natural distribution to obtain this 
kind of sampling is setting the density of sub-bands proportional 
to a power law of frequency, V , where a is a given (negative) 
constant, i.e, 



r y > i - 1 f VM 



Therefore, 



y i = \(<-< 1 ) 



/- i 

N-l 



+ v 



,d+i 



l/ta+l) 



(10) 



For the special case of a - - 1 we have instead 



v, = exp (In v M - In v m ) 



i-\ 

N-l 



+ In v„ 



(11) 



Equation [TUl becomes Eq.|8]for a = 0. It can be shown that 
the standard deviation of vj l (i.e., the precision in the estimate of 
the ionospheric dispersion) is maximum for a = -5 /3 in the case 
of a large number of sub-bands. We show in Fig.|2]the optimum 
a as a function of the number of sub-bands, N. These values of 
a were computed by finding numerically the minimum of cr(K) 
(using Eq. |6]l as a function of a. The values of a shown in Fig. 
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Fig. 2. Boxes, optimum values of a (i.e., those that minimize the 
formal uncertainty in the estimate of the ionospheric dispersion) 
as a function of the number of sub-bands. Line, model corre- 
sponding to Eq.rT2l 

[2] can be estimated using the model (which is also shown in the 
figure) 



0.6 



a = 



JV-4.35 



5 
3' 



(12) 



which tends to -5/3 for large N. These are the values of a that 
were used to obtain the results reported in the next section. We 
show an example of this power-law distribution in Fig.QJc). 

3.4. Golomb rulers 

A Golomb ruler is a set of n, integers such that the set of differ- 
ences, dij - rij-rij, has no repeated elements (see, e.g., Atkinson, 
Santoro & Urrutia 1986). It is intuitive that Golomb rulers are a 
good choice to maximize the dynamic range of the fringes, since 
all pairs of sub-bands are separated incoherently one respect to 
the other. Therefore, the sidelobes of the Fourier transform of the 
bandpass are minimum. The improvement in the fringe dynamic 
range when the sub-bands are distributed according to Golomb 
rulers has been previously reported for the case of 8 sub-bands 
(Mioduszewski & Kogan 2004). Here we generalize the study 
to different number of sub-bands and also analyze the impact of 
this kind of distribution in the precision of the estimate of the 
ionospheric dispersion. The central frequencies, v,, of the sub- 
bands are computed using the equation 



1+/3— y„ 

n N 



(13) 



where n,- is the zth element of a Golomb ruler of N elements (by 
convention, n\ - 0). The Golomb rulers for N < 24 were taken 
from the OGR project at http://distributed.net/ogr, and the others 
from Atkinson, Santoro & Urrutia ( 1986). We show an example 
of this distribution in Fig.QJd). 



obtain information from as many regions of the bandwidth as 
possible. Moreover, a more homogeneous sampling makes eas- 
ier the the connection of the phases between the sub-bands, since 
an unsampled wide lag in the spectrum could contain a num- 
ber of 2n phase cycles that could introduce biases in the data 
analysis. To better understand this statement, let us consider, for 
example, a non-dispersive delay, r, added to the fringe. The dif- 
ferential phase between sub-bands i and j, due to t, would be 

A<f> u = 4>j - (pi = 2ttt(vj - v,) 

For larger values of Vj - v, (i.e., for wider lags between 
sub-bands), the differential phase between sub-bands is larger. 
Therefore, the probability of the differential phase to be larger 
than 2tt is higher for wider lags between sub-bands. 

From this point of view, the uniform and/or the power-law 
distributions would be the best frequency configurations for the 
interferometer. However, we must also take into account the dy- 
namic range of the fringes and the quality in the estimate of the 
atmospheric dispersion, which are analyzed in the following sub- 
sections. 



4.1. Ionospheric dispersion 

In Fig.[3fa), we show the uncertainty in the estimate of the iono- 
spheric dispersion, cr(K), in units of the minimum possible un- 
certainty (i.e., o-(K) m j„, computed from Eqs.|6]and|7]l, as a func- 
tion of the number of sub-bands, for a total bandwidth of [5 — 1.0 
(see Eq. Q]). The four different distributions are shown. For all 
distributions, the uncertainty (relative to the minimum possible 
one) increases with the number of sub-bands. It can be seen that 
the distribution based on Golomb rulers give the most precise es- 
timates of the ionospheric dispersion, followed by the power-law 
distribution (with an uncertainty ~ 3% larger, depending on the 
number of sub-bands), and the random and uniform distributions 
(with uncertainties ~ 6% larger, also depending on the number 
of sub-bands). 

In Fig.|31 we show the uncertainty in the estimate of the iono- 
spheric dispersion, in units of the minimum possible uncertainty, 
as a function of the bandwidth (i.e., /?, see Eq. []]) using a total 
of 32 sub-bands to cover the bandwidth. Golomb rulers yield 
again the most precise estimates of the ionospheric dispersion, 
although the uncertainty increases as the bandwidth increases. 
On the contrary, the power-law distribution keeps the uncertainty 
roughly constant as a function of the bandwidth. Both the con- 
stant and random distributions also increase the uncertainty in 
the atmospheric dispersion as the bandwidth increases, being 
this uncertainty ~5% larger than that obtained with the Golomb 
rulers. 

We conclude that the power-law distribution and that based 
on Golomb rulers give higher precisions in the estimate of the 
ionospheric dispersion. Although the difference between uncer- 
tainties from all the distributions is not so large (lower than 10% 
in all cases), its optimization may be important to obtain high- 
contrast images. 



4. Results 

Figure[T]shows that the random and Golomb-rulers distributions 
tend to poorly sample some regions of the spectrum and over- 
sample others. On the contrary, the constant and the power-law 
distributions sample the bandwidth in a more homogeneous way. 
A more homogeneous sampling of the spectrum is preferable to 



4.2. Fringe dynamic range and delay ambiguity 

In Fig.[3jb), we show the dynamic range of the fringes as a func- 
tion of the number of sub-bands for a total bandwidth of f3 — 1 .0 
(see Eq. [TJ. We notice, however, that the dynamic range of the 
fringes is independent of/3, since a change in/3 is equivalent to a 
change in the delay scaling of the fringes (regardless of a phase 
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Fig. 3. (a) Formal uncertainty in the estimate of the ionospheric 
dispersion, in units of the minimum possible uncertainty, (b) dy- 
namic range of the fringes, and (c) delay ambiguity (i.e., dis- 
tance between the fringe peak and the closest sidelobe) in units 
of the Nyquist time resolution. All these quantities are shown as 
a function of the number of sub-bands for /? = 1 .0 (see Eq. [TJ. 

factor that depends on v m ). Two regions in the sub-band space 
can be readily seen. 

For N < 20, the random and the power-law distributions give 
higher dynamic ranges. Suprisingly, for these values of N, the 
distribution based on Golomb rulers give dynamic ranges ~ 1 . 
Why? Golomb rulers are sets of integer numbers. Therefore, the 
Fourier transforms of these sub-band distributions are periodic. 
If the delay window is larger than the period of the Fourier trans- 
form, there will be more than one peak in the fringe. The fringe 
period depends on each ruler and increases with the number of 
channels. For N < 20, the fringe period is shorter than our delay 
window (1024 times the Nyquist time resolution, i.e., 512 chan- 
nels in each direction of the delay) so there is more than one peak 
in the fringe. For N > 20, the fringe period is larger and only 
one peak remains in the delay window. For N > 20, Golomb 
rulers give the highest dynamic ranges (around 20 - 30% higher 
than those based on the random and power-law distributions). 



1.8 

jl.75 
b" 1.7 



« 



1.65 



1.6 




0.5 



a 



1.5 



Fig. 4. Formal uncertainty in the estimate of the ionospheric dis- 
persion, in units of the minimum possible uncertainty, as a func- 
tion of /3 (see Eq.[T) for the case of 32 sub-bands. 



In all cases, the uniform distribution gives very poor dynamic 
ranges, ~ 1, as it is indeed expected, since the Fourier transform 
of the bandpass is a periodic function with a very short period. 
For the case of the delay ambiguity, strong changes are seen as 
a function of the number of sub-bands for the Golomb rulers 
(the ambiguity ranges between ~ 20 and 512 channels) and the 
power-law distribution (the ambiguity ranges between ~ 150 and 
512 channels, although the lower limit increases with N). These 
changes in the delay ambiguity are due to several sidelobes with 
similar peak values. Changing N also changes the relative height 
of the sidelobe peaks. As a consequence, for different values of 
N, different sidelobes are selected as the highest and, therefore, 
very different delay ambiguities are obtained. On the contrary, 
the random distribution has a delay ambiguity of ~ 250 channels 
for all values of N (we notice, however, that, for this distribution, 
the figure shows the average of 100 different fringes). The uni- 
form distribution, as expected, has a very small delay ambiguity 
(lower than 100 channels). This ambiguity increases with N, also 
as expected, since the spacing of sub-bands (which is shorter for 
larger N) is inversely proportional to the period of the fringe. 

A first conclusion is that the uniform distribution is not a 
good choice from the point of view of the quality in the esti- 
mate of the group delay. The Golomb rulers are a good choice 
when the number of sub-bands is large (N > 20, although this 
number decreases if the width of the delay window decreases). 
The power-law distribution is, in general, a good choice for 
all N. It gives the best compromise between homogeneity in 
the sampling of the bandwidth, precision in the estimate of the 
ionsopheric dispersive effects, dynamic range of the correlated 
fringes, and ambiguity of the group delay. Therefore, this would 
be the preferable sub-band distribution to use in low-frequency 
(wide-band) interferometric observations. 

Nevertheless, these conclusions are based on a number of 
sub-bands up to 64. If the number of sub-bands is large (say, 
N = 512) there are no available Golomb rulers to work with, 
but we can still compare the results obtained from the uniform, 
random, and power-law distributions. 

Setting N = 512, the dynamic range of the fringes is simi- 
lar for the three distributions, if we use a delay window of 1024 
Nyquist channels (using 512 sub-bands, the period of the fringe 
corresponding to the uniform distribution is longer than the de- 
lay window). However, the power-law distribution still gives 
lower formal uncertainties in the estimate of the ionospheric dis- 
persion (around 10% lower than the other distributions for/3 = 2 
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Fig. 5. Formal uncertainty in the estimate of the source spectral 
index, y, in units of that of the uniform distribution (<x(-y) uni ), as 
a function of the bandwidth, /J (see Eq. Q}, for the case of 64 
sub-bands. 



and 4% for/? =1). Therefore, the power-law distribution is still 
the best choice with a number of sub-bands as large as 512. 



4.3. Source spectral index 

Wide-band observations allow to precisely determine the spec- 
tral indices and spectral curvatures of radio sources. The dif- 
ferent distributions of sub-bands may also affect the achievable 
precision in the estimate of the spectral properties of the radio 
sources. For the case of the spectral index, y (being the flux 
density, S oc v y ), the formal uncertainty, cr(-y) depends on the 
distribution v,- in the form 



0"(y) = 0" (log v,) 



(14) 



We show in Fig.[5]the formal uncertainty in the estimate of y, 
in units of that of the uniform distribution, cr(y) uni , as a function 
of the bandwidth, /J, for the case of 64 sub-bands. It can be seen 
in the figure that the uniform, random, and power-law distribu- 
tions give similar precisions in the estimate of the spectral index 
(although the precision for the power-law distribution slightly 
increases with the bandwidth). Surprisingly, the Golomb-ruler 
distribution gives a precision ~ 5% higher than the other distri- 
butions for all bandwidths. 

The formal uncertainties of the uniform, random, and power- 
law distributions for a large number of sub-bands {N = 512) 
are similar to those shown in Fig. [5] Therefore, for wide-band 
observations using many sub-bands (i.e. where no Golomb rulers 
are available), the use of the power-law distribution is the best 
choice, allowing for a 1 - 2% higher precision in the estimate 
of the spectral index of the sources. With a smaller number of 
sub-bands, the use of Golomb-ruler distributions would improve 
the precision in the estimate of y by ~ 5%. 

4.4. Other contributions to the optimum power-law 
distribution 

Other contributions to the visibility phases (either due to the 
electron plasma of the ionosphere or to chromaticity in the struc- 
ture of the observed sources), as well as the contribution of the 
galactic radiation to the visibility noise, have not been consid- 
ered in the previous sections. In this section, we study how 
these contributions may affect the optimum sampling of the 



ionospheric dispersion using the power-law distribution of sub- 
bands. 



4.4.1 . Plasma frequency of ionospheric electrons 

Equation [5] holds in the region of frequencies much higher 
than the plasma frequency, v p , of the ionospheric electrons. The 
electron density in the ionosphere takes values around 10 4 - 
10 6 e~ cm 3 . This translates into a plasma frequency in the range 
~ 1 - 10 MHz (e.g., Pacholczyk 1970 Eq. 2/72). From the re- 
fraction index of a plasma (e.g., Pacholczyk 1970, Eqs. 2.77 and 
2.78), the phase drift in the case of v, ~ v p is 



A(pi - K'vi 



V,(V; + V B ) 



- 1 



(15) 



where vb is the Larmor electron frequency (for the Earth mag- 
netic field it takes the value ~ 1 MHz) and the ± signs corre- 
spond to the two possible circular polarizations of the radiation. 
We notice that the effect of vb in the computations reported in 
this section is negligible. K' is related to the TEC above each 
element of the baseline in the line-of-sight to the source. 

We computed the optimum values of a (i.e., the exponent 
of the power-law distribution of sub-bands) that optimize the 
sampling of the ionospheric phase drifts (i.e., minimize the for- 
mal uncertainty in the estimate of K') for observing frequencies 
close to v p . We call these values a p , to distinguish them from 
the values without the effect of the plasma frequency (i.e., a, 
which are shown in Fig. [2] and given in Eq. fT2l . In Fig. [6] we 
show the ratios a p /a as a function of the minimum observing 
frequency, v m , in units of the plasma frequency, v p . For instance, 
using a bandwidth of /J = 1.0 (i.e., vm = 2v m ) a plasma fre- 
quency v p = 10 MHz, and a minimum frequency v m = 100 MHz 
(i.e, 10 times the plasma frequency) results in values of a p equal 
to 0.86 times those of a. For a large number of sub-bands (i.e., 
for a = -5/3) this results in a p = -1.43. 

We notice that the exponent a p approaches zero as the mini- 
mum frequency approaches to the plasma frequency. Therefore, 
the power-law distribution approaches the constant distribution 
at very low frequencies. Increasing the bandwidth tends to com- 
pensate a bit the decrease in the absolute value of a p (i.e., the 
ratios a p /a increase when /3 increases) although this effect is 
small, as it can be appreciated in Fig. [6] 

4.4.2. Frequency-dependent source structure 

A source structure which is independent of the observing fre- 
quency is a strong assumption over the broad frequency ranges 
considered in the previous sections. The contribution of a possi- 
ble source chromaticity in the visibility phases can be divided in 
two parts. One is related to the source being intrinsically differ- 
ent at different frequencies (this contribution might be especially 
important for extended sources) and the other one is related to 
the position of the source (or that of its brightest feature) being 
different at different frequencies (as it is the case, for instance, 
of a self-absorbed core-jet structure). 

On the one hand, the contribution of the source structure to 
the visibility phases can be determined from the image of the 
source at each frequency. This image depends on the visibil- 
ity calibration, but can also be used to refine such calibration. 
Therefore, it should be possible, in principle, to decouple the 
source structure (which introduces baseline-dependent phases) 
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Fig. 6. Optimum values of a p , in units of a, as a function of the minimum observing frequency, v m , in units of the plasma frequency, 
v p . In each figure, the values are computed for five different values of v p (2, 4, 6, 8, and 10 MHz). Each figure corresponds to a 
different bandwidth, /3 (0.1, 0.5, 1.0, and 1.5). 



from the ionospheric dispersion (which introduces antenna- 
dependent phases), with the help of iterative and elaborated 
imaging-calibration algorithms. The details of these algorithms 
and their impact in the precision of the estimated ionospheric 
contribution, as it is decoupled from the source-structure contri- 
bution after the imaging, is beyond the scope of this paper. 

On the other hand, the contribution of a frequency-dependent 
source position on the visibility phases can be studied if some as- 
sumptions are considered. Porcas (2009) reported on the effects 
of source chromaticity in the source position estimates through 
VLBI astrometry, performed using either phase delays or group 
delays, for the case of a core-jet structure following the model of 
Blandford & Konigl dl9791 l. The contribution of the chromatic 
core-shift to the interferometric phase is 



A0 str = K s v 1 



(16) 



where K s depends on the physical conditions in the jet and the 
angle of the projected baseline with respect to that of the jet. 
The parameter 6 also depends on the physical conditions in the 
jet and may take values between and ~ 2. If we combine 
Eqs. 131 and [TBI we obtain the total drift of the visibility phases 
through the sub-bands when both effects, ionospheric dispersion 
and source chromaticity, are taken into account: 



A(f>i = (K + K s v^ s )v' 1 . 



(17) 



Therefore, the effect of a frequency-dependent source po- 
sition is equivalent to the addition of an extra term coupled to 
the parameter K of the ionospheric dispersion. For 6 = (i.e., 
no core-shift), K s translates into a contribution to (f> equal to a 



constant (i.e., non-dispersive) group delay, with the phase pro- 
portional to the frequency. This group delay is equivalent to a 
shift of the source, which is the same at all frequencies. This 
shift can be easily fixed in the calibration if the source position 
is known. For 6-1, the effect of the source chromaticity on <p 
is just adding a constant, so it does not affect cr(K) (see Eq. [6J. 
Therefore, the sub-band distribution for the optimum estimate of 
ionospheric dispersion will be the same as with no chromatic ef- 
fect. However, for the (physically unrealistic) case 6 - 2, there is 
a complete coupling between the ionospheric dispersion, K, and 
the source chromaticity, K s , so it is especially difficult to cal- 
ibrate the ionospheric delay using these sources, regardless of 
the distribution of sub-bands used. Real sources may have val- 
ues of 6 falling between the values here analyzed (0, 1, and 2). It 
is thus expected to find intermediate cases in which the effect of 
the frequency-dependent source position will be a combination 
of either a non-dispersive group delay, no effect in the estimate 
of ionospheric dispersion, and a complete coupling between that 
estimate and the source position. Additionally, all the core-jet 
sources detected in a given FoV may have different values of 
6, so different calibration issues will appear in the same image 
depending on the coordinates of each source and the values of 6. 
In any case, we notice that K s depends on the direction of the 
projected baseline relative to that of the jet, so it is a baseline- 
dependent quantity. However, the ionospheric contribution, K, 
is antenna-dependent. This different behavior of K s and K, de- 
pending on the pair of stations selected, should allow for a robust 
decoupling of K s from K, provided the number of observing sta- 
tions is large enough. Therefore, any chromaticity in the source 
structure and/or position should not affect the results reported in 
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this paper, as long as the baseline-dependent chromatic effects 
are decoupled from the antenna-dependent ionospheric contri- 
bution using the appropriate calibration algorithms. 

4.4.3. Radiation from the Galaxy 

For frequencies below ~ 400 MHz, the sky brightness temper- 
ature is dominated by the Galactic radiation, which depends 
strongly on the observing frequency (T^ y oc v y with y ~ 2.5). It 
means that in LOFAR wide-band observations, the noise in the 
low-frequency sub-bands will be higher than that in the high- 
frequency sub-bands. In the cases of observations dominated by 
the radiation from the Galaxy, Eq.|6]must be adapted to take into 
account the different uncertainties in each sub-band. 

Thermal noise from the sky brightness temperature translates 
into a gaussian-like noise in the real and imaginary parts of the 
visibilities, with a value of <x proportional to the equivalent flux 
density of the system, which is in turn proportional to the total 
(i.e., receivers plus source) temperature (e.g., Thomson, Moran, 
& Swenson |1991| ). If the observed sources are strong, the noise 
in the amplitudes and phases can also be approximated as be- 
ing gaussian-like, with a <x proportional to that of the real and 
imaginary parts of the visibilities. If we take this approximation 
into account and assume that the galactic radiation dominates 
the system temperature (i.e., r sys ~ r s k y ), then the uncertainty 
in the visibility phase of the zth sub-band is cr, oc vr 5 and Eq. [6] 
becomes 



0.5 



cr(K) 



^ 



EvfZvf-(Evf) 2 



(18) 



The values of a in Eq. [10] that minimize the cr(K) given by 
Eq.[l8](we call these values a ga i) are shown in Fig. [7] as a func- 
tion of the normalized bandwidth, /3. For small values of /?, we 
find that a ga \ is positive. It means that the noise in the phases 
at the lower frequencies is such high, that a better sampling of 
the high-frequency region of the band yields more precise esti- 
mates of the ionospheric dispersion. However, this effect is less 
important for wider bands and/or higher minimum frequencies, 
as it can be seen in Fig. [7] If the band is wide enough (for in- 
stance, yS > 1.7 for a minimum frequency of 200 MHz) the opti- 
mum value ag a i is negative. The 20 - 90 MHz band of LOFAR 
corresponds to a value of j3 as large as 3.5, for which we find 
a gal = -0.8 (setting v m = 20 MHz). 

5. Summary 

We studied the achievable precision in the estimate of the iono- 
spheric dispersion, the ambiguity of the group delay, the dy- 
namic range of the correlated fringes, and the precision in the 
estimate of the source spectral index in low-frequency and wide- 
band interferometric observations for four different distributions 
of the sub-bands through the total bandwidth of the detectors: 
constant spacing between sub-bands, random spacings, spac- 
ings based on a power-law distribution, and spacings based on 
Golomb rulers. 

For a large number of sub-bands, spacings based on Golomb 
rulers give the most precise estimates of dispersive effects and 
the highest dynamic ranges of the fringes. Spacings based on 
the power-law distribution give similar (but slightly worse) re- 
sults, although the results are better than those from the Golomb 
rulers for a smaller numbers of sub-bands. Random distributions 
of sub-bands result in relatively large dynamic ranges of the 
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Fig. 7. Optimum value of a that minimizes the cr(K) given in 
Eq.[l8]for different values of the minimum frequency, v m . These 
values of a g . d \ have been computed assuming that the Galaxy ra- 
diation dominates the noise of the visibilities. A line at a g - d \ = 
is also shown for clarity. 



fringes, but the estimate of dispersive effects through the band 
is worse. A constant spacing of the sub-bands results in very bad 
dynamic ranges of the fringes, but the estimates of dispersive ef- 
fects have a precision similar to that obtained with the power-law 
distribution. 

From all combinations of the number of sub-bands and the 
total covered bandwidth, the power-law distribution (with a 
given by Eq. fT2l gives the best compromise between homo- 
geneity in the sampling of the bandwidth, precision in the esti- 
mate of the ionsopheric dispersive effects, dynamic range of the 
correlated fringes, and ambiguity of the group delay. Therefore, 
this would be the preferable sub-band distribution to use in low- 
frequency (wide-band) interferometric observations. 

Finally, we study how the power-law distribution that opti- 
mally samples the ionospheric dispersion is affected in the cases 
of 1) observing frequencies close to the plasma frequency of the 
ionospheric electrons, 2 chromatic effects in the structure of the 
sources, and 3) non-negligible noise coming from the Galaxy 
radiation. 
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